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Abstract 

Background: Sick leave due to neck pain (NP-SL) is costly and negatively impacts the productivity of the nursing 
and midwifery workforce. Identification of modifiable risk indicators is necessary to inform preventive efforts. This 
study aimed to investigate the role of pain-related psychological features (pain catastrophizing, fear of movement, 
and pain coping) in NP-SL alongside other potential risk indicators. 

Methods: A cross-sectional analysis of a large cohort study of Australian and New Zealand nurses and midwives, 
established between 1st April 2006 to 30th March 2008, was undertaken. Recruitment procedures adopted within 
each Nursing Council jurisdiction were governed by the individual regulatory authorities and their willingness to 
engage with the study. Invitations directed potential participants to a purpose-built internet-based survey, 
where study information was provided and consent requested. Once consent was obtained, a range of standardized 
tools combined into one comprehensive electronic questionnaire was elicited. Exposure variables assessed included 
pain characteristics and a broad range of psychological, psychosocial, occupational, general health and demographic 
factors. Two-way interactions between age and gender and candidate exposures were also assessed. Binary logistic 
regression was performed using manual backward stepwise elimination of non-significant terms. 

Results: The cohort included 4,903 currently working nurses or midwives aged 18-65 years. Of these, 2,481 
(50.6%) reported neck pain in the preceding 12 months. Our sample comprised of 1,854 working nurses and 
midwives with neck pain in the preceding year who supplied sick leave data. Of these, 343 (18.5%) reported 
taking sick leave in the preceding year due to their neck pain. The final most parsimonious multivariable model 
demonstrated neck pain severity (adjusted odds ratio, [aOR] = 1.59), passive pain coping (aOR= 1.08) and fear of 
movement (aOR = 1 .06) increased the likelihood of NP-SL in the previous year. Interactions between demographic 
and general health factors exhibited both protective and risk relationships with NP-SL, and there was no association 
between pain catastrophizing and NP-SL. 

Conclusions: Findings demonstrate that sick leave due to neck pain was associated with pain severity, fear of 
movement and passive pain coping. In addition, there were complex interactions found between demographic 
and general health factors. These features represent potentially modifiable targets for preventive programs. 
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Background 

Neck pain (NP) appears to be more prominent in the 
general population than previously considered [1]. It af- 
fects up to 71% of individuals during their lifetime [2], 
and general population studies indicate that between 30- 
50% of individuals will experience NP in any 12-month 
period (range: 12.1%-71.5%) [3]. Furthermore, it is com- 
monly recurrent, with over 50% of workers with NP 
reporting symptoms a year later [4]. Economic analyses 
reveal that NP can place a considerable annual financial 
burden on society [3,5,6]. In the Netherlands, sick leave 
due to neck pain (NP-SL) was reported to comprise over 
a quarter of the total NP financial burden, and costed an 
estimated US$185.4 million from 1.4 million days of lost 
productivity associated with worker absence [5]. 

Given the societal, workforce and personal implica- 
tions of NP-SL, the factors associated with NP-SL must 
be understood in order to inform preventive efforts. Ex- 
posures examined to date have predominated around 
pain characteristics and occupational, demographic and 
psychosocial factors with inconsistent results [7-13]. In 
addition, there has been only one study of factors associ- 
ated with NP-SL as an exclusive outcome [13]. This ana- 
lysis showed that work-related neck flexion and neck 
rotation, low decision authority and medium skill discre- 
tion were significantly related to NP-SL in industrial 
workers [13]. All other studies have examined work ab- 
senteeism due to pain in composite regions that in- 
cluded the neck in combination with the shoulder, upper 
extremity and/or upper and lower back [7-12]. Given 
that factors related to composite outcomes may differ to 
those related to NP-SL in isolation, analyses of sick leave 
due to pain in the neck region alone are needed to clar- 
ify exposure relationships. 

The relationship between NP-SL and pain-related psy- 
chological features has never been investigated. Theoretical 
models posit that mal-adaptive pain-related cognitions and 
behaviours contribute to development of chronic pain 
[14,15]. Specifically, fear- avoidance theory asserts that cata- 
strophizing about pain leads to pain-related fear (e.g. fear 
of movement) and avoidance/escape behaviours (e.g. pas- 
sive pain coping) that lead to disability, disuse, and further 
pain [14,15]. Staying home from work is a potential proxy 
of avoidant behaviour, hence mal-adaptive pain-related fea- 
tures may play a role in explaining NP-SL. The demon- 
strated relationship between pain-related psychological 
features and sick leave due to low back pain [16,17] sug- 
gests these factors may be similarly related to NP-SL. 
Given that psychological features may be modifiable, and 
therefore ideal targets for preventive efforts, their role in 
NP-SL must be defined. 

Utilizing a large occupational cohort of Australian and 
New Zealand nurses and midwives, the aim of this study 
was to determine the role of pain-related psychological 



factors (pain catastrophizing, fear of movement and pain 
coping) and other psychosocial, occupational, demo- 
graphic and pain characteristics in explaining NP-SL. 
Under the Australian National Employment Standards, 
full-time employees are entitled to 10 days' paid personal 
leave (for sick and paid carer's leave) per year. Part-time 
employees receive a pro rata entitlement to sick leave 
based on the number of hours they work. In New Zealand 
there is a minimum provision of five days' paid sick 
leave a year after the first six months of continuous em- 
ployment and an additional five days' sick leave after 
each subsequent 12-month period. Leave accumulates 
from year to year in both countries. The nursing popu- 
lation was targeted to answer our research questions 
since NP negatively impacts this profession. A popula- 
tion study of multiple occupational groups in Britain 
found that nurses reported the highest annual preva- 
lence (15%) for NP that prevented normal activities 
[18]. Further, workforce attrition studies show than 
nurses with neck/shoulder pain are 1.5 times more 
likely to leave nursing work [19]. Hence, an understand- 
ing of the factors associated with NP-SL can inform pro- 
fessional policy and preventive programs to maximise 
the viable nursing workforce. 

Methods 

Participants, study design, and procedure 

The sample comprised Australian and New Zealand 
nurses and midwives recruited as part of the Nurses and 
Midwives e-Cohort Study (NMeS). The NMeS is a longi- 
tudinal study that investigates workforce and health 
outcomes in nurses and midwives from Australia, New 
Zealand and the United Kingdom (UK) using innova- 
tive online (http://www.e-cohort.net) data collection 
methods [20]. Recruitment strategies, response rates 
and characteristics of the NMeS sample have previ- 
ously been described in detail [21,22]. In brief, invitations 
for participation varied between the varying State-based 
and national Nursing and Midwifery Council jurisdictions, 
and included: flyers within licence renewal packages 
(Queensland, Tasmania, New Zealand), postcards mailed to 
home addresses (Australian Capital Territory, New South 
Wales, Western Australia), and mass advertisement and/or 
articles available and circulated to Council members 
(Northern Territory, South Australia, Victoria). Invitations 
directed potential participants to a purpose-built internet- 
based survey (see: http://www.e-cohort.net). Upon entering 
the study web-site, detailed information was provided and 
consent requested. Once received, participants established 
their personal profile and completed the registration page 
which requested baseline demographic, contact and work- 
related details. Participants were then directed to the 
baseline survey, which contained various standardized 
instruments around the Work/Life Balance and Staying 
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Healthy study themes and took approximately 40 minutes 
to complete. Participants with incomplete surveys were 
emailed reminders approximately two weeks, four weeks 
and six weeks following initial registration within the 
NMeS. Recruitment, registration and the baseline survey 
was open from 1 April 2006 to 31 March 2008. Follow-up 
surveys have taken place approximately bi-annually: the 
first opened on 30 August 2008 and closed 26 September 
2009, and the second opened on 25 May 2010 and closed 
on the 31 May 2011. The cross-sectional data used in this 
analysis was taken from the NMeS baseline survey col- 
lected in 2006-2007. Participants included in this analysis 
were aged 18-65 years, currently employed in nursing or 
midwifery work, had professional registration in Australia 
or New Zealand and reported NP in the preceding 
12 month period. As UK nurses and midwives were re- 
cruited using opportunistic strategies, as opposed to the 
targeted approach used in Australia and New Zealand, 
they were omitted from this study to minimize responder 
bias. Overall, 7,812 (2.3%) eligible nurses and midwives 
participated from a total pool of 334,400. However, age, 
gender, occupational and health profiles were similar be- 
tween countries and to national workforce figures [21,22]. 

Ethical approval for the NMeS was granted clearance 
by the University of Queensland's Behavioural and Social 
Science Ethics Review Committee (No. 2005000696) 
and the Massey University Human Ethics Committee 
(Wellington: No. 05/71), and all participants provided 
informed consent. 

Outcome measure 

The outcome for all analyses was whether NP-SL was 
taken in the previous year, using a binary (Yes/No) vari- 
able. A question from the General Standardised Nordic 
Questionnaire was used to determine 12-month NP 
prevalence [23]. Neck trouble was defined as 'ache, pain 
or discomfort in the area extending from the base of the 
skull to the top of the trunk and shoulder region'. 

Exposure measures 

Pain, psychosocial, demographic, occupational, and gen- 
eral health factors were all investigated. 

Pain characteristics 

Neck pain severity (rating of worst pain on an 11-point 
scale, 0-10) was measured with an item from the 
Graded Chronic Pain Scale adapted to inquire about the 
preceding 12 month interval [24]. 

Psychosocial factors 

Four psychosocial domains identified by Linton [25] that 
are relevant to neck and back pain were assessed: cogni- 
tive, emotional, social and behavioural. 



(i) Cognitive: Pain catastrophizing is an exaggerated 
negative orientation to pain [26] that includes 
"performing cognitive activities that exacerbate the 
fearful aspects of the pain experience" [27]. We 
used the reliable and valid Pain Catastrophizing 
Scale to assess this construct, using its three 
subscales rumination (range 0-16 points), 
magnification (range 0-12 points), and helplessness 
(range 0-24 points) [26,28]. Fear of movement was 
assessed using the 11-item version of the Tampa 
Scale of Kinesiophobia (TSK-11) (range 11-44 
points) that has demonstrated reliability and 
validity in patients with chronic low back pain 
[29] . The scale assesses the construct kinesiophobia, 
which is defined as "an excessive, irrational and 
debilitating fear of physical movement and activity 
resulting from a feeling of vulnerability to painful 
injury or re-injury" [29]. Respondents were asked 
to answer TSK-11 questions specific to their NP 
experience. We also asked participants whether 
they believed nursing work caused or exacerbated 
their NP, and quantified this belief with a binary 
(Yes/No) variable we named Belief that nursing 
work caused/exacerbated NP. 

(ii) Emotional: The 10-item version of the Centre 
for Epidemiologic Studies Depression Scale 
(CESD-10) was used to assess depression [30]. A 
current depressive symptoms binary (Yes/No) 
variable was generated using the cut-off score 

of >10 [30]. Prior diagnosis of anxiety and diagnosis 
of depression by a medical professional were 
assessed with binary (Yes/No) variables. A 
general mental health status measure was 
included using the Mental Health Component 
Score from the SF-36 (Version 2) Health Survey, 
standardised to 1998 United States general 
population normative values [31]. Higher scores 
on this scale reflect better health. 

(iii) Social: The Job Content Questionnaire (JCQ) 
(Version 1.1) was used to assess the following 
work-related constructs: job control (the sum of 
Skill Discretion + Decision Authority subscales), 
social support (the sum of Co-Worker Support + 
Supervisor Support subscales) and psychological 
demands [32]. The Effort-Reward Imbalance 
Questionnaire (ERIQ) was used to assess effort- 
reward imbalance (ERI) and over-commitment 
(OC) [33]. Possibly due to challenges of electronic 
administration of the ERIQ [22], 69% of ERI data 
was missing and consequently ERI was not included 
in analyses. OC measures the inability to withdraw 
from work obligations [33], and is categorized into 
high (upper tertile of the sample distribution) and 
low (lower two tertiles) [34] . 
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(iv) Behavioural: pain coping was assessed with the 11-item 
abbreviated version of the Vanderbilt Pain Management 
Inventory (VPMI) [35,36] that has established internal 
consistency and factorial validity in individuals with 
spinal pain [37]. The instrument generates two 
subscales: active coping and passive coping. The 
active coping subscale assesses such strategies as 
functioning in spite of pain or distracting oneself 
from pain (score 5-25 points) and the passive coping 
subscale assesses behaviours such as depending on 
others for pain management assistance and allowing 
pain to adversely impact usual activities (score 6-30 
points). As for the TSK-11, respondents were 
requested to answer VPMI questions specific to 
their NP experience. 

Demographic factors 

The demographic factors included in the model were 
gender, marital status, country and age, with age stan- 
dardised to years at 30th November 2006. 

Occupational factors 

A categorical variable named job classification was used to 
quantify nursing and midwifery work into five groups that 
delineate qualifications and clinical/ non-clinical work: (i) 
registered nurses and midwives; (ii) enrolled nurses and 
nurse assistants; (iii) specialist registered nurses and mid- 
wives (e.g. nurse practitioners, clinical educators); (iv) 
administrators, managers and policy makers; and, (v) 
academics and university educators. Current shift work 
was assessed with a binary variable (Yes/No). Manual 
handling task frequency was assessed with an 8-item 
specially-developed instrument (total score 0-24 points). 
The instrument inquires about frequency (assessed with a 
Likert scale 0, 1-4, 5-10, 11+ that is correspondingly 
scored 0-3) of manual handling activities such as assisting 
patients to transfer from sit to stand and pushing/pulling 
lifting equipment. General physical demands was assessed 
with the JCQ subscale of the same name that includes two 
questions relating to physical effort and rapid and con- 
tinuous physical activity. 

General health factors 

A range of general health factors that could potentially 
have a confounding influence were included in analyses. 
Current pregnancy, hours of sleep and smoking status 
variables were included, and the Physical Health Com- 
ponent Score of the SF-36v2, standardised to 1998 US 
population normative data [31], was used as a measure 
of general physical health. Finally, a range of diagnosed 
conditions (asthma, pneumonia, osteoarthritis, rheuma- 
toid arthritis, chronic obstructive pulmonary disease, 
multiple sclerosis, Parkinson's disease, stroke, breast 
cancer, congestive heart failure, high blood pressure, 



lung cancer) were considered in analyses using binary 
(Yes/No) variables. 

Test-retest reliability 

Measures used within the NMeS baseline survey were 
sourced from existing instruments with established psy- 
chometric properties when available. When question- 
naires were unavailable or inappropriate, questions were 
constructed or adapted and assessed for test-retest reli- 
ability in a sub-study of 32 nurses and midwives. These 
nurses and midwives re-completed the baseline NMeS 
survey on average 8.8 ± 2.5 days (range 5-16) after their 
first completion. The reliability of categorical data was 
assessed with proportion of observed agreement (p a ) 
and the kappa statistic (k). Specifically, the maximum 
kappa obtainable (K max ) was calculated, and the propor- 
tion of K max achieved (K/K ma x) was used to represent 
test-retest reliability [38]. The reliability of continuous 
variables was assessed with the intra-class correlation 
coefficient using two-way models with random effects 
and single measures (ICCp, ij) [39]. The reliability of 
variables assessed in this way was acceptable: annual sick 
leave due to neck pain, p Q - 0.81, i</K max = 0.74; belief 
that nursing work caused/exacerbated NP, p a = 0.82, 
K/i<m a x = 0.39; pain severity ICC[2, i] = 0.73 (95% confi- 
dence interval [CI]: 0.42, 0.89); and manual handling 
task frequency ICC [2 , 1} = 0.98 (95% CI: 0.95, 0.99). 

Statistical analysis 

The characterisation of some candidate exposure vari- 
ables could be expressed in multiple forms (e.g. continu- 
ous, categorical, dichotomous, as well as individual 
versus composite variables). Selection between these 
characterisations were made in bivariable logistic regres- 
sion models using the Bayesian Information Criterion 
(BIC) [40]. The BIC penalises for model complexity and 
rewards for goodness-of-fit; with the preferred model 
balancing these competing demands and yielding the 
lowest BIC value. On this basis, for example, the pain 
catastrophising subscales were included in multivariable 
modeling over the total pain catastrophising score, and 
the composite JCQ variables (i.e. job support and job 
control) were selected rather than individual variables. 
Crude relationships between the outcome and all expos- 
ure variables were then assessed using bivariable logistic 
regression. Since confounding can lead exposures with 
modest bivariable associations to have significant rela- 
tionships in multivariable models [41], a conservative 
threshold for inclusion in multivariable analysis was set 
at p < 0.40. Next, all exposure variables yielding a bivariable 
association with p < 0.40 were included in a multivariable 
main effects logistic regression model. Non-significant ex- 
posure variables were then sequentially omitted from the 
model if they failed to reach then p < 0.40 criterion. This 
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process determined the most parsimonious main effects 
model. At this stage, age and gender were retained within 
the statistical model, regardless of significance. Next, all 
two-factor interaction terms between age and gender and 
all remaining exposure variables were entered into the 
model. Stepwise elimination of non-significant interaction 
terms was then undertaken using a p < 0.05 criterion. Non- 
significant main effect variables that did not form part of a 
significant interaction term were then also eliminated, in a 
stepwise fashion, using the p < 0.05 until the final model 
was ascertained. The goodness-of-fit of the final model was 
assessed with the Hosmer-Lemeshow test [42]. Multicolli- 
nearity between regressors was evaluated by tolerance and 
condition index estimates derived from linear regression 
models. All analyses were undertaken using Stata version 
12.0 (StataCorp, College Station, Texas, USA). 

Results 

Sample characteristics 

The NMeS recorded 10,721 participant registrations, 
which equated to 10,112 unique registrants following re- 
moval of duplicates (some individuals registered on 
more than one occasion). Of 8,651 individuals who went 
on to complete the online survey instrument, 7,812 were 
professionally registered in Australia or New Zealand 
and 4,903 were aged 18-65 years and were currently 
working in a nursing or midwifery role. Of these, 50.6% 
of these individuals reported NP in the preceding 
12 months. Our study sample for the present analysis is 
comprised of 1,854 individuals with NP in the preceding 
year that provided data regarding sick leave (i.e. NP-SL) 
during that period. Compared to those who responded 
to the NP-SL question, non-responders were older 
(45.0 years verses 43.9 years, p = 0.004), more likely to 
be female (95.4% versus 92.9%, p = 0.03), and more 
likely to be Australian (82.4% versus 88.8%, p < 0.001), 
but no different in marital status (p = 0.35) or job 
classification (p = 0.08) profiles. The differential com- 
pletion of the sick leave response question by eligible 
participants who reported NP in the preceding 12 months 
may introduce bias in pursuant prevalence estimates and 
reported associations. 

Sick leave due to neck pain characteristics 

Overall, 343 (18.5%) nurses and midwives our study 
sample reported taking sick leave in the preceding year 
due to their NP, while 1,511 (81.5%) took no leave due 
to this pain. Table 1 presents the demographic, occupa- 
tional and general health characteristics of the sample 
partitioned by the dichotomous sick leave status, to- 
gether with the bivariable associations between these 
factors and NP-SL. Based on the p < 0.40 exclusion cri- 
terion, variables measuring co-morbid low back pain, 
pregnancy, and smoking status were not considered 



within the multivariable model development. Variables 
measuring age and gender were retained at this stage. 

Similarly, the pain and psychosocial characteristics of 
the sample together with the bivariable associations are 
presented in Table 2. The only factor that did not meet 
the inclusion threshold for the multivariable analysis was 
psychological demands. All remaining variables were 
considered in the multivariable model development. 

Multivariable analysis 

In developing the most parsimonious main effects logis- 
tic regression model, the sequentially eliminated vari- 
ables (in order) were: active coping; manual handling 
task frequency; country; over commitment; hours of 
sleep; diagnosed depression; asthma; belief that nursing 
work caused/exacerbated NP; high blood pressure; mari- 
tal status; job classification; physical demands; and de- 
pressive symptoms. Introducing the two-factor (age and 
gender) interaction terms and eliminating non-significant 
variables in a stepwise fashion yielded the adjusted odds 
ratio (aOR) estimates, together with the associated 95% 
CIs, presented in Table 3. 

There were significant main effect relationships be- 
tween NP-SL and neck pain severity, passive coping, and 
fear of movement. However, more complicated signifi- 
cant relationships emerged between NP-SL and age and 
general physical health scores, and between NP-SL and 
gender and general mental health scores; see Table 3. 
There was no reason to reject the adequacy of this final 
model (Hosmer-Lemeshow test, p = 0.96), and the model 
explained 31% of the NP-SL variance. 

Post hoc analysis 

A post hoc analysis was performed to examine why pain 
catastrophizing was not retained in the final model. Pain 
catastrophizing was moderately correlated with passive 
coping (R = 0.39), and had modest correlations with the 
other significant variables: pain severity (R = 0.12); gender 
{R = -0.01); age (R = -0.11); general physical health score 
{R = -0.24); general mental health score (R = -0.25); and 
fear of movement {R = 0.24). Linear regression analysis of 
pain catastrophizing yielded tolerance estimates that 
ranged from 0.63 (passive coping) to 0.98 (age x general 
physical health score interaction) and a condition index 
estimate of 13.6, suggesting that multicollinearity among 
independent variable in not a significant issue. 

In the bivariable analysis with NP-SL, pain catastro- 
phizing was significant (p < 0.001) but explained only 2% 
of the variability of NP-SL. When a trivariable analysis 
was performed with passive coping added to the model, 
pain catastrophizing was no longer statistically signifi- 
cant (p = 0.26), and the model explained 14% of NP-SL 
variability. The interaction of pain catastrophizing and 
passive coping on NP-SL yielded an estimated coefficient 
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Table 1 Demographic, occupational and general health characteristics of Australian and New Zealand nurses and 
midwives with neck pain in the preceding 12-months by whether they took no sick leave (n = 1511) or sick leave (n = 343), 
together with the bivariable associations between these factors and neck pain related sick leave (NP-SL) 

No sick leave taken Sick leave taken OR (95% CI) p-value 

Demographic factors 



Age in years, mean (SD) 


44.3 


(9.3) 


42.1 


(9.8) 


0.98 


(0.96, 0.99) 


<0.001 


Gender, n {%) 














0.52 


Female 


1407 


(93.1) 


316 


(92.1) 


1.00 


(reference) 




Male 


104 


(6.9) 


27 


(7.9) 


1.16 


(0.74, 1 .80) 




Marital status, n (%) 














0.37 


Never married 


239 


(16.1) 


64 


(18.9) 


1.00 


(reference) 




Married/de facto 


1049 


(70.4) 


235 


(69.3) 


0.84 


(0.61, 1.14) 




Divorced/separated/widowed 


201 


(13.5) 


40 


(11.8) 


0.74 


(0.48, 1.15) 




Country, n (%) 














0.01 


New Zealand 


283 


(18.7) 


44 


(12.8) 


1.00 


(reference) 




Australia 


1228 


(81.3) 


299 


(87.2) 


1.57 


(1.11, 2.20) 




Occupational factors 
















Job classification, n (%) 














0.36 


Registered nurse/midwife 


700 


(67.0) 


143 


(60.6) 


1.00 


(reference) 




Enrolled nurse/nurse assistant 


85 


(8.1) 


27 


(114) 


1.56 


(0.97, 2.49) 




Specialist registered nurse/midwife 


125 


(12.0) 


32 


(13.6) 


1.25 


(0.82, 1 .92) 




Administrator/manager 


98 


(9.4) 


25 


(10.6) 


1.25 


(0.78, 2.01) 




Academic/university lecturer 


36 


(3.4) 


9 


(3.8) 


1.22 


(0.58, 2.60) 




Current shift work, n (%)* 


825 


(54.8) 


214 


(62.8) 


1.39 


(1.09, 1.77) 


0.01 


Manual handling frequency, mean (SD) 


5.7 


(6.0) 


6.9 


(7.0) 


1.03 


(1.01, 1.05) 


0.002 


Physical demands, mean (SD) 


4.9 


(1.6) 


5.1 


(1.6) 


1.08 


(1.00, 1.16) 


0.04 


General health factors 
















General physical health, mean (SD) 


50.9 


(7.2) 


47.5 


(8.3) 


0.95 


(0.93, 0.96) 


<0.001 


Co-morbid low back pain, n (%)* 


1061 


(71.5) 


237 


(70.4) 


1.06 


(0.82, 1 .37) 


0.68 


Pregnant, n (%)* 


35 


(2.4) 


6 


(1.8) 


0.76 


(0.32, 1 .83) 


0.54 


Asthma, n (%)* 


327 


(22.9) 


95 


(29.5) 


1.41 


(1.08, 1.85) 


0.01 


Emphysema/chronic bronchitis, n (%)* 


30 


(2.1) 


3 


(0.9) 


0.44 


(0.13, 1.50) 


0.18 


Pneumonia, n (%)* 


141 


(9.8) 


25 


(7.8) 


0.77 


(0.50, 1 .20) 


0.25 


High blood pressure, n (%)* 


221 


(15.4) 


62 


(19.3) 


1.31 


(0.96, 1 .79) 


0.09 


Hours of sleep, n (%) 














0.06 


> 8 


270 


(26.1) 


52 


(22.3) 


1.00 


(reference) 




6-7 


675 


(65.3) 


150 


(64.4) 


1.15 


(0.82, 1 .63) 




< 5 


88 


(8.5) 


31 


(13.3) 


1.83 


(1.10, 3.03) 




Smoking status, n (%) 














0.93 


Never smoked 


749 


(51.9) 


120 


(51.8) 


1.00 


(reference) 




Ex-smoker 


489 


(33.9) 


108 


(33.1) 


0.98 


(0.75, 1 .28) 




Current smoker 


206 


(14.3) 


49 


(15.0) 


1.05 


(0.74, 1 .50) 




Note: "denotes categorical variables which have an omitted reference group of 'no'. 



of 0.0002 (Standard error [SE]: 0.0019, Wald's \ statistic Discussion 

p = 0.91). These analyses suggest that the bivariable The key objective of this study was to determine the dif- 

effect of pain catastrophizing on NP-SL was primarily ferential role of pain-related cognitions and behaviours 

confounded by passive coping. alongside other risk indicators in NP-SL. While many 
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Table 2 Pain and psychosocial characteristics of Australian and New Zealand nurses and midwives with neck pain in 
the preceding 12-months by whether they took no sick leave (n = 1511) or sick leave (n = 343), together with the 
bivariable associations between these factors and neck pain related sick leave (NP-SL) 

No sick leave taken Sick leave taken OR (95% CI) p-value 

Pain characteristics 



Neck pain severity, mean (SD) 


5.2 


(2.2) 


7.5 


(1.8) 


1.70 


(1.59, 1.83) 


<0.001 


Psychosocial factors: cognitive 
















Fear of movement, mean (SD) 


20.8 


(5.0) 


24.2 


(5.1) 


1.15 


(1.12, 1.18) 


<0.001 


Pain catastrophising, mean (SD) 


6.4 


(7.4) 


8.4 


(8.3) 


1.03 


(1.02, 1.05) 


<0.001 


Rumination subscale, mean (SD) 


2.6 


(3.0) 


3.0 


(3.2) 


1.04 


(1.00, 1.08) 


0.03 


Magnification subscale, mean (SD) 


1.3 


(1.7) 


1.6 


(1.8) 


1.11 


(1.04, 1.19) 


0.001 


Helplessness subscale, mean (SD) 


2.5 


(3.4) 


3.7 


(4.0) 


1.08 


(1.05, 1.12) 


<0.001 


Belief work caused/exacerbated NP, n (%)* 


1232 


(82.1) 


312 


(91.5) 


2.35 


(1.57, 3.51) 


<0.001 


Psychosocial factors: behavioural 
















Pain coping - passive coping, mean (SD) 


1 1 .4 


(3.6) 


14.2 


(4.0) 


1.21 


(1.17, 1.25) 


<0.001 


Pain coping - active coping, mean (SD) 


16.7 


(3.4) 


16.5 


(3.1) 


0.98 


(0.95, 1 .02) 


0.33 


Psychosocial factors: emotional 
















Depressive symptoms, n (%)* 


347 


(24.2) 


118 


(37.2) 


1.86 


(1.43, 2.40) 


<0.001 


Diagnosed anxiety, n (%)* 


232 


(16.2) 


73 


(22.7) 


1.52 


(1.13, 2.05) 


0.01 


Diagnosed depression, n (%)* 


393 


(27.4) 


119 


(37.1) 


1.56 


(1.21, 2.02) 


0.001 


General mental health, mean (SD) 


46.1 


(10.1) 


42.8 


(10.6) 


0.97 


(0.96, 0.98) 


<0.001 


Psychosocial factors: social 
















Skill discretion, mean (SD) 


37.6 


(4.7) 


37.1 


(5.1) 


0.98 


(0.96, 1 .00) 


0.10 


Decision authority, mean (SD) 


36.0 


(6.9) 


35.4 


(7.4) 


0.99 


(0.97, 1.01) 


0.21 


Job control, mean (SD) 


73.6 


(10.4) 


72.6 


(11.4) 


0.99 


(0.98, 1 .00) 


0.12 


Psychological demands, mean (SD) 


16.0 


(3.1) 


16.1 


(3.2) 


1.01 


(0.98, 1 .04) 


0.62 


Co-worker support, mean (SD) 


6.0 


(1.0) 


5.9 


(1.0) 


0.92 


(0.81, 1.03) 


0.15 


Supervisor support, mean (SD) 


5.8 


(1.5) 


5.6 


(1.7) 


0.90 


(0.84, 0.98) 


0.01 


Job support, mean (SD) 


11.9 


(2.2) 


11.5 


(2.4) 


0.94 


(0.89, 0.99) 


0.02 


Over commitment, n (%) 














0.004 


Low 


729 


(69.9) 


142 


(60.2) 


1.00 


(reference) 




High 


314 


(30.1) 


94 


(39.8) 


1.54 


(1.15, 2.06) 





Note: *denotes categorical variables which have an omitted reference group of 'no'. 



cases of sick leave due to incapacitating NP are likely to 
be appropriately exercised, fear-avoidance theory [15] 
suggests there may also be cases where fearful cognitions 
about pain prohibit work attendance rather than actual 
functional limitation. It is this NP-SL that has the poten- 
tial for prevention if relevant psychological factors are 
identified and targeted interventions developed. In a co- 
hort of nurses and midwives with NP in the preceding 
year, we found pain severity, fear of movement and pas- 
sive pain coping increased the likelihood of work ab- 
sence. In addition, there were complex interactions 
found between demographic and general health factors. 

The relationship between neck pain severity and NP- 
SL was expected. Pain severity has consistently shown a 
relationship with NP-SL [43-45], and in our model, 



severity of worst pain in the preceding year increased 
the odds of NP-SL by 59% per unit increase in pain rat- 
ing. Age and gender have also been significantly associ- 
ated with NP-SL [45,46], although this study is the first 
to demonstrate significant effect modifications with gen- 
eral physical and mental health scores. This may par- 
tially explain some of the equivocal findings reported in 
earlier studies [46]. 

We found that fear of movement increased the odds 
of NP-SL in the preceding year by 6% per unit increase 
in TSK-11 score, suggesting that even after controlling 
for pain severity, nurses and midwives with NP who are 
fearful of movement are more likely to stay home from 
work. This is an exciting finding given that fear of move- 
ment has shown some potential for modification in an 
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Table 3 Adjusted odds ratio (aOR) estimates and 
associated 95% CIs for the terms in the final 
multivariable logistic regression model 



Variables 


aOR 


(95% CI) 


Age (in years) 


0.975 


(0.959, 


0.990) 


General physical health score 2 


0.996 


(0.975, 


1.017) 


Age 1 x General physical health score 2 


1.002 


(1.000, 


1 .004) 


Neck pain severity score 


1.593 


(1470, 


1.727) 


Passive coping score 


1.079 


(1.033, 


1.126) 


Fear of movement score 


1.057 


(1.022, 


1.093) 


General mental health score 2 


0.994 


(0.979, 


1 .009) 


Male 3 


0.876 


(0.453, 


1.695) 


General mental health score 2 x Male 3 


0.937 


(0.890, 


0.987) 



1 Centred around 45 years (median value); 2 centred around 50 (population 
norm score); 3 taking females as the reference category. 



intervention study with individuals following back injury 
[47]. Further, the current finding is consistent with evi- 
dence from other studies that demonstrate fear of 
movement is related to sick leave due to low back pain 
(LBP-SL) [16,17]. In a previous analysis of LBP-SL in 
nurses and midwives, fear of movement increased the 
likelihood of LBP-SL by 5% for women and 8% for men 
[16]. Hence, the strength of association in this study is 
consistent with previous analyses for LBP-SL. In terms 
of Hill's criteria for causation [48], there is mounting 
evidence for a causal relationship between fear of 
movement and sick leave due to spinal pain: the temporal 
relationship [17], strength and consistency of association 
in this and others studies [16,17], and plausibility of the 
relationship [15]. However, before causation can be 
asserted, all of Hill's criteria need to be satisfied, in- 
cluding the yet to be established dose-response rela- 
tionship between fear of movement and sick leave due 
to spinal pain, demonstrated prevention or amelior- 
ation by appropriate experimental regimen, and the es- 
tablishment of coherency and adequate specificity. 

Passive coping also increased the likelihood of NP-SL, 
by 8% per unit increase in score (on a 24-point scale). 
This finding adds weight to the importance of passive 
coping in musculoskeletal pain outcomes. Mercado and 
colleagues [49] demonstrated that individuals using mod- 
erate to high levels of passive coping strategies were five 
times more likely to develop disabling neck and/or low 
back pain than those with low levels. Passive coping has 
also been shown to increase the likelihood of LBP-SL by 
7% per unit increase in score [16]. 

The significant interaction terms are illuminating. The 
first reveals that as age and general physical health 
scores together increase, the odds of NP-SL decrease less 
than if age and general physical health scores were con- 
sidered separately as independent effects. Similarly, as 
age and general physical health scores together decrease, 



then the odds of NP-SL increase less than if considered 
separately. So, in workforce planning, both variables 
should be considered together rather than in isolation. 
The second interaction term reveals that the general 
mental health score association with NP-SL is primarily 
associated with males; and the lower the general mental 
health score a male has, the higher his odds of NP-SL. 
While significant here, interaction terms are rarely 
investigated. 

Our analysis identified no direct relationship between 
pain catastrophizing and NP-SL. Our a priori hypothesis 
that pain catastrophising would be associated with NP- 
SL was based on fear-avoidance theory [14,15] which as- 
serts pain catastrophizing leads to fear of movement and 
avoidant behaviours. While significant in the bivariable 
model, NP-SL was primarily confounded by passive cop- 
ing in the multivariable model. This finding, and a similar 
non-significant relationship between pain catastrophising 
and LBP-SL [16], suggests that pain catastrophizing may 
be related to proximal factors in the fear- avoidance model 
but not directly related to work absence. Another rejected 
a priori hypothesis was that of an association between ac- 
tive coping and NP-SL. We found active coping had no re- 
lationship with NP-SL in either bivariable or multivariable 
analyses. A lack of association has also been demonstrated 
between active coping and NP-related disability [16,49]. 
According to fear-avoidance theory, active coping re- 
sponses should represent confrontational behaviours that 
lead to recovery [14,15]. Given that flaws in the measure- 
ment of active coping have been proposed [49], further 
work is needed to establish valid and reliable measure- 
ment of active coping behaviours to facilitate further in- 
vestigations into the possible protective influence of active 
coping on NP-SL. 

Self-reported psychosocial work exposure such as job 
strain, low co-worker support, decreased job security 
and overall stress at work have also been shown to be 
risk factors for NP [45,46]. While many considered oc- 
cupational and social variables had significant bivariable 
associations with NP-SL in this study, none appeared in 
the final multivariable model. This implies, for our study 
group, that demographic, general health, pain characteris- 
tics, and psychosocial factors mediated these occupational 
and social variables in explaining NP-SL. The relatively 
homogeneous and structured occupational and social 
characteristics around nursing and midwifery may hide 
the small but significant associations seen elsewhere [46] . 

An urgent call has been made for high quality rando- 
mised controlled trials of workplace interventions for 
neck pain [45]. Our findings advocate for multidimen- 
sional preventive interventions that incorporate strat- 
egies to ameliorate both fear of movement and passive 
coping behaviours. Fortunately, a coping-oriented cogni- 
tive and behavioural intervention developed by Linton 
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and colleagues [50,51], is clearly defined and could read- 
ily be applied in future trials in nurses and midwives. 
However, that intervention did not reduce fear of move- 
ment [50,51]. Program strategies that target fear of 
movement, such as those included in the Pain-Disability 
Prevention Program [47], could be combined with the 
coping-oriented program in the hope of additive effects 
on NP-SL. Optimized pain management should be also 
considered in multidimensional NP-SL preventive inter- 
ventions since pain severity is related to NP-SL and has 
also been shown to predict future fear of movement 
[52]. Given that pain severity, fear of movement and pas- 
sive coping are related to both NP-SL (as describe here) 
and sick leave due to low back pain [16], a multidimen- 
sional program to address these factors is likely to have 
positive impacts on work loss due to both neck and low 
back pain. 

This study has a number of strengths but also limita- 
tions. This is the first investigation of pain-related psy- 
chological factors and NP-SL and provides new insights 
into exposure relationships. Pain coping has been con- 
sidered as a confounder in an analysis of sick leave due 
to a composite outcome of neck or upper extremity pain 
[12], but the main effect of pain-related psychological 
characteristics has never before been explored. There 
are three important limitations to consider: (i) the cross- 
sectional design precludes causal inferences; (ii) all NMeS 
data were self-reported and do not include objective mea- 
sures; and, (iii) the NMeS sample represents less than 5% 
of the nursing and midwifery population that investigators 
attempted to invite [21,22]. 

Only with the following of NMeS participants over 
time, together with the acquisition of empirical informa- 
tion gleaned from other studies, will causal inferences 
able to be made. Self-report studies have validity prob- 
lems; participants may exaggerate conditions or events 
in order to make their situation seem worse or through 
social desirability bias they might under-report condi- 
tions or events in order to conceal or minimize stigma- 
tised matters. Participants may also simply be mistaken 
or misremember the material covered by the survey. 
However, our deliberate choice of standardised instru- 
ments sought to minimise this inherent effect. 

The modest participation rate potentially threatens the 
external validity of the data, although the associated im- 
pact of differential participation in large-scale, multiple- 
domain, targeted population studies is likely to be small 
[22]. Selective participation in relation to NP is unlikely 
because potential participants were informed in NMeS 
promotional material about broad research objectives 
(i.e. to study workforce and health outcomes) that 
made no reference to NP. Moreover, in relation to 
basic demographic characteristics and NP prevalence, 
there is no obvious responder bias [21,22]. In this 



analysis, 50.6% of nurses and midwives who met inclu- 
sion criteria reported NP in the previous year. This is 
reassuringly consistent with a sample of New Zealand 
nurses with similar demographic characteristics who 
reported 52% prevalence using the same measurement 
instrument [53]. Further, in our sample 18.5% of nurses 
with NP reported NP-SL in the preceding year, compared 
to 21.3% of individuals in a large general population sam- 
ple that reported neck, shoulder or upper back pain [54]. 

The overall NMeS registrants have been shown to be 
similar to the source Australian and New Zealand popu- 
lations in relation to gender, age and Indigenous repre- 
sentation [21,22]. This similarity is notable as it might 
have been opined that older participants would have re- 
duced access to, and have been less comfortable with, a 
fully electronic survey. However, the use of information 
technology is now considered a core aspect of modern 
practice and is a basic element in nursing and midwifery 
programmes [20]. Furthermore, Internet access at home 
is becoming commonplace. In 2008, 75% of Australian 
households had a computer and 67% had home Internet 
access, of which 78% was broadband [20]. Internet pene- 
tration is likely higher in our target populations as rates 
increase markedly with higher incomes. In respect to the 
sample of nurses and midwives included in this analysis, 
participants' demographic characteristics (92.9% female 
and mean 43.9 years) were consistent with the Austra- 
lian nursing and midwifery labour force at the time of 
recruitment (90.4% female and mean 43.7 years) [55]. 
Slight differences may reflect demographic associations 
with NP prevalence, since NP appears to differentially 
affect women (females report NP more than men in 83% 
of prevalence studies) [2]. 

Conclusions 

In summary, pain severity, fear of movement and passive 
coping behaviours are associated with sick leave due to 
NP in nurses and midwives. Other factors that relate to 
NP-SL include demographic and general health factors. 
Interestingly, occupational, social, and several cognitive 
factors (including pain catastrophizing) were not signifi- 
cantly associated with NP-SL in our sample. Our find- 
ings indicate that multidimensional interventions to 
reduce fearful cognitions about pain, ameliorate passive 
coping behaviours and optimise pain management 
should be developed and examined via high quality trials. 
Further research will determine whether such interven- 
tions can reduce the societal, workforce and individual 
costs of NP-SL. 
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